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Spherical collapse of the Bose-Einstein Condensate (BEC) dark matter model is studied in the 
Thomas Fermi approximation. The evolution of the overdensity of the collapsed region and its 
expansion rate are calculated for two scenarios. We consider the case of a sharp phase transition 
(which happens when the critical temperature is reached) from the normal dark matter state to the 
condensate one and the case of a smooth first order phase transition where there is a continuous 
conversion of “normal” dark matter to the BEC phase. We present numerical results for the physics 
of the collapse for a wide range of the model’s space parameter, i.e. the mass of the scalar particle 
m x and the scattering length l s . We show the dependence of the transition redshift on m x and l s . 
Since small scales collapse earlier and eventually before the BEC phase transition the evolution of 
collapsing halos in this limit is indeed the same in both the CDM and the BEC models. Differences 
are expected to appear only on the largest astrophysical scales. However, we argue that the BEC 
model is almost indistinguishable from the usual dark matter scenario concerning the evolution of 
nonlinear perturbations above typical clusters scales, i.e., > 10 14 Mq. This provides an analytical 
confirmation for recent results from cosmological numerical simulations [H.-Y. Schive et al, Nature 
Physics, 10 , 496 (2014)]. 
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I. INTRODUCTION 

It is widely accepted that dark matter is one of the 
main components of the universe. Due to the strong 
observational evidence corroborating its existence many 
different areas of physics have incorporated dark mat¬ 
ter related investigations in their agenda. According to 
the standard cosmological model, dark matter composes 
around 1/4 of the universe’s energy budget and 5/6 of the 
total matter. Baryons represent the remaining fraction of 
the latter. This picture has been confirmed by different 
data, but remarkably by the latest Planck results jl|. 

The crucial aspects of these studies concern the parti¬ 
cle nature and the astrophysical/cosmological behavior of 
such component. At the particle level, candidates belong¬ 
ing to the WIMP (weakly interacting massive particles) 
category produce a viable model (see [§]—and references 
therein -for a very recent review). Also, for the homo¬ 
geneous, isotropic and expanding background, the dark 
matter ensemble should present a vanishing pressure in 
order to enable structure formation Q. 

Although the success of the Cold Dark Matter (CDM) 
scenario it is important to mention some of its draw¬ 
backs. The theoretical clustering patterns (calculated 
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via numerical simulations) of CDM particles at galactic 
level correspond to the NFW profile [IJ which is cuspy 
at the centre of the particle distribution. This seems 
to be in clear contradiction to the observed velocities 
in the central region of galaxies which demand a cored 
distribution. At the same time, the simulated distribu¬ 
tion of satellites around typical Milk Way like galaxies 
shows one order of magnitude excess of sub-structures 
which are not observed. These two issues are known as 
the cusp-core problem and the missing satellites problem, 
respectively. Even if baryonic physics in such simulations 
could eventually alleviate these problems, it is not clear 
so far whether or not CDM is the correct model for the 
dark matter phenomena. See []5] and references therein. 

One can argue that dark matter is a pathological man¬ 
ifestation of choosing Einstein’s general relativity (GR) 
as the gravitational theory. This suspect is the pillar of 
a research line in which modified gravity theories are in¬ 
voked. See jdj for reviews on modified gravity models and 
their observational constraints. However, reliable exper¬ 
iments at the solar system level confirm GR predictions 
with great accuracy [7]. Therefore, this fact seems to 
be powerful enough to keep in a fist moment GR as our 
standard description for gravitational interaction. 

Since there are no confirmed evidence to abandon GR, 
dark matter remains being essential and therefore one 
needs new alternatives within this context. In this case, 
the possibilities are also vast. The classical ones were hot 
dark matter (HDM) [§] and warm dark matter (WDM) 
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0. While the former has been ruled out due to the pos¬ 
itive observation of galaxies below the jeans mass scale 
of relativistic dark matter particles, the latter is one of 
the leading rivals of CDM. Indeed, particles with masses 
m ~ keV fit the WDM spirit. They are not as light as 
HDM particles and therefore allowing the existence of 
structures and, at the same time, not as heavy as CDM, 
in such a way that there would exist some suppression 
mechanism able to alleviate the small scale problems of 
the CDM paradigm (see however |lQ] for a recent discus¬ 
sion of WDM results). Models with a similar clustering 
dynamics as WDM are, for instance, fuzzy dark matter 
M , the self interacting dark matter [l 2| and the viscous 
dark matter 0. 

In this work we study a dark matter model which has a 
different nature. Let us assume 0—spin DM particles hav¬ 
ing therefore a bosonic distribution. As predicted and al¬ 
ready observed in laboratory bosonic particles are able to 
condensate 0 (see also jjj|), occupying the same energy 
state and forming the so-called Bose-Einstein conden¬ 
sates when their temperature reaches the critical value 
T crt . Of course, this phenomenon occurs under very con¬ 
trolled experimental situations, but one might wonder 
in principle what happens if the same would happen on 
astrophysical scales. 

Although quite hypothetical this description could 
serve as an effective approach for understanding dark 
matter as a cosmological scalar field </> whose dynamics 
is driven by some repulsive potential V(<p). This gives 
rise to the Bose-Einstein Condensate (BEC]_ dark matter 
model which has been widely studied |16H2ll| . The main 
idea is that normal, i.e., non-condensate, dark matter un¬ 
dergoes a phase transition at some critical redshift z cr t 
during the universe’s evolution. Then, independently on 
the details of the transition, all the dark matter converts 
into the condensate state forming a BEC “fluid” 1 * . 

The dynamics of BEC systems is studied via the Gross- 
Pitaevskii equation, which is a nonlinear Schrodinger 
equation [23|. From this starting point, the Madelung 
decomposition is used to transform the BEC dynamics 
into a set of fluid equations resulting in an effective pos¬ 
itive pressure. With such fluid picture one is able to in¬ 
vestigate astrophysical/cosmological problems. This pro¬ 
cedure will be shown in more details in the next section. 

The general aspects of this model concerning the back¬ 
ground evolution and the linear perturbations are already 
very well understood [2^ - (30i] . But, in order to fully un¬ 
derstand the final clustering patterns of the BEC dark 
matter model high resolution hydrodynamical/N-body 
simulations are still needed & More recently, Ref. 
0 has formulated smoothed-particle hydrodynamics 
numerical methods to solvinggeneral Gross-Pitaevskii- 
Poisson system. Schive et al [0] provided recently high- 
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resolution cosmological simulations for the model. They 
obtained that there is a remarkable difference at the in¬ 
ternal galactic level, i.e., its density profile. The latter 
result is indeed desired. However, they found that BEC 
DM is indistinguishable from CDM at large cosmological 
scales. Our focus here in this work is to understand such 
latter claim. From the theoretical point of view, a first 
step on this issue is the study of the nonlinear gravita¬ 
tional collapse in a cosmological background. Concerning 
the BEC dark matter model, recently Ref. [Hj addressed 
the collapse of “already formed BEC condensates”, i.e., 
only the post-transition stage. Nevertheless, a realistic 
configuration can be more complicated since it also in¬ 
volves the dynamics of the baryonic component as the 
universe evolves from the matter to the dark energy dom¬ 
ination epochs. Moreover, the phase transition can also 
take place during the evolution of the collapsed region. 
Therefore, especially for galaxy cluster scales, the evolu¬ 
tion of the background cosmological dynamics should be 
taken into account. 

We will perform in this work a natural extension of 
Ref. 0] which has analysed the “free-fall” collapsed 
of a BEC dark matter sphere. However, we assume a 
more realistic cosmological scenario where dark matter 
coexist with baryons and a cosmological constant. Then, 
we address the correct case where the transition occurs 
during the nonlinear clustering process. 

Fundamental quantities here are the condensate pa¬ 
rameters, namely, the mass of the particle m x and the 
scattering length l s . They determine the moment at 
which the phase transition takes place z cr t and the speed 
of sound in the condensate fluid, for example. After the 
critical redshift z cr t one can admit two different dynam¬ 
ics. The simplest case is to assume an abrupt transition, 
i.e., for z < z cr t all dark matter obeys the Bose-Einstein 
dynamics. This seems to be a reasonable approximation 
to the problem. This situation will be studied in section 

uni 

One can also assume the case in which the full con¬ 
version of all dark matter occurs in a finite time and it 
finishes at a redshift Zbec < z cr t■ Therefore, the phase 
transition lasts a finite time in which a mixture of “nor¬ 
mal” and condensate dark matter make up the total mat¬ 
ter component. We study in section HVl this case. 

We present our results covering many order of magni¬ 
tude in the model parameter space 10 -6 meV< m x < 10 4 
meV;10 -12 fm< l s < 10 12 fm. Interesting quantities to 
be found here are the final (at 2 = 0) value of the density 
contrast and the expansion rate and the redshift of the 
turnaround z a , i.e., the moment at which the collapsed 
region detaches from the background. 

In summary, this paper has the following structure. In 
the next section we develop the background dynamics of 
the BEC dark matter. We present in section Hill general 
equations for the spherical top-hat collapse formalism. 
These equations will be studied in more detail in sec¬ 
tions m and 0 where, respectively, we address the case 
of abrupt transition and the usual phase transition. We 
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conclude in the final section. 


II. THE BACKGROUND DYNAMICS OF THE 
BOSE-EINSTEIN CONDENSATE DARK 
MATTER 


In this work we always have a flat background dynam¬ 
ics composed by baryons, dark matter and a cosmological 
constant. This expansion rate reads 

H 2 = — 2 — ( Pb + Pdm + Pa) • (1) 

The post-decoupling dynamics of the baryonic compo¬ 
nent is assumed to be pressureless Pb = 0 and there¬ 
fore pb = Pbo(l + z ) 3 where pbo is its density today 
at z = 0. Its value is such that pbo = U-boPcO, where 
p c o = 3Hq/8ttG. We can safely adopt f^o = 0.05 accord¬ 
ing to nucleosynthesis constraints. The Hubble constant 
assumed here is Hq = 70 Km/s/Mpc. We will also fix 
fldmo = 0.25 or equivalently Ha = 0.75. 

The difference here from the standard ACDM model 
will be the dark matter dynamics. Before the transi¬ 
tion takes place, at temperatures T > T crt ; or redshifts 
z > z cr t, DM behaves as an isotropic gas in thermal 
equilibrium. From kinetic theory the pressure of a non- 
relativistic gas in this regime is given by 

9s f q 2 c 2 . 3 g s [ q 4 o 

w " = 3F, / ~ md p “ ^sF J m 

( 2 ) 

with a 2 = (v 2 ) /3c 2 , where g s is the number of spin de¬ 
grees of freedom, h the Planck constant, q the momentum 
of a particle with energy E = and distri¬ 

bution function /. A typical value for the velocity disper¬ 
sion is a = 3 x 10 -6 . In practice, since this quantity can 
be seen as the dark matter equation of state parameter 
u>dm = Pdm/Pdm this value is consistent with the assump¬ 
tion of pressureless fluid usually adopted for CDM. Note 
that the full relativistic fluid is obtained when (u 2 ) = c 2 . 

After dark matter’s conversion it obeys the conden¬ 
sate dynamics, which is governed by the Gross-Pitaevskii 
equation 

H 2 o 

ih— = -—V 2 * + V(r,t)'b+g( (3) 

at 2to x 

where m x is the mass of the particle and V (r, t) is the 
trapping potential. The non-linearity term with only 
two-body interparticle interaction (quadratic) reads 

g(m (4) 

where Uo = 4irh 2 l s /m x . This definition has the funda¬ 
mental parameters of the model, namely the scattering 
length l s and the particle mass m x . The former is asso¬ 
ciated to the nature of the short range self-interactions 
in the condensate. For example, in laboratory systems, 
it can be either positive (the case of Rb 87 atoms with 


l s = 5.45 nm and then repulsive interactions) jUJ or 
negative (the case of Li 7 atoms with l s = 1.45 nm and 
then attractive interactions) [36]. In this work we will 
consider only cases where / s > 0. The impact of l s on 
the mass-radius configurations of astrophysical BEC has 
been investigated in [371 ]. 

Note that there appears some degeneracy for the Uq 
parameter, i.e., there are infinities combinations of l s and 
m capable to produce the same Uo value. We discuss this 
degeneracy and the admissible numerical values of these 
parameters in the next sections. 

In order to apply the Gross-Pitaevskii equation to 
astrophysical problems one proceeds with the so-called 
Madelung decomposition. In this procedure, the wave 
function is replace by 


$ = V®ei s ( r ’ ( > (5) 

where p = | V I'| 2 is the number density of the system and 
S is the velocity potential. The mass/energy density can 
be written in terms of the mass of each individual particle 
as p x = m x p. 

Therefore, the BEC system can be described in terms 
of a hydrodynamical set of equations, which are 


du 




_ +{u-S7)u= -V 

at p x \ m 


V 


9p x 


+ v ■ (Px^) = 0 > 


dt ' ' 
where we define 


VQ 


m 


. ( 6 ) 


(7) 


u = -VS, 


m y 


Q = - 


h 2 

2m x VK 


( 8 ) 

(9) 


The particle self-interaction of this specific BEC- 
inspired fluid gives rise to a pressure of polytropic form 


Px = 


27 rh 2 l s 2 

^ P\' 
m x K 


( 10 ) 


On the other hand, the quantum potential Q/m x results 
in the often called quantum pressure 2 . We can use the 
identity 

djPij _ (%Q (-q) 

Px m x 

where is the quantum anisotropic pressure tensor [37j , 
given by 


h 2 (dipdjp 

P * = M l — 



( 12 ) 


2 Note that both quantum pressure and self-interaction pressure 

are of quantum mechanical origin. 
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For the problem we have in mind, the potential V ( r , t) 
in © is in fact the gravitational potential which is 
sourced by p x via the Poisson equation V 2 F = 47 t Gp x . 
This allows us to solve the system of equations. 

In the cases where the pressure due to the particle self¬ 
interaction dominates, the quantum anisotropic pressure 
can be neglected. This is the so called Thomas-Fermi ap¬ 
proximation. In |38| the authors estimated, for the case 
of BEC dark matter halos, in which cases the Thomas- 
Fermi limit is valid. They consider the forces associated 
with both pressures, that balances the gravitational col¬ 
lapse, and find that the Thomas-Fermi regime is valid 
when 

— >2, (13) 

where 

k = Airh 2 -^— . (14) 

m x 

Adopting R as the mean radius and M as the mass of a 
BEC dark matter halo it is found that 


2 h2 R 

Kh = —TTn - . 

3 M 


(15) 


The quantity (1151) can be written in characteristic val¬ 
ues 


k h = 2.252 x 10 


-64 


R 


100 kpc 


10 12 


M 


— ) eV cm 3 . 


(16) 

If we consider halos with a size between the Milky Way 
(M = 10 12 M q and R = 100 kpc) and a typical dwarf 
galaxy (M = 10 10 Mq and R = 10 kpc) we can constraint 
kh in the range 


k h ^ 2 x (10“ 64 - 10~ 63 ) eV cm 3 . (17) 


A cosmological dark matter fluid with the above pres¬ 
sure leads to the background expansion 

H 2 = ~(Pb + Px + Pa) ■ (19) 

where p x is the BEC dark matter density, which in the 
Thomas-Fermi limit is determined by the pressure USD 
via the continuity equation. 

More details will be discussed in sections HVlandlVl In 
fact, we will follow in this work the background expansion 
determined in ref. [l8i]. 


III. THE NONLINEAR TOP-HAT COLLAPSE 

Here we present the basic equations that describe the 
evolution of a spherical collapsing matter region in an 
expanding background. This is the ideal technique for 
studying the clustering patterns of dark matter halos. 

We will follow standard calculations presented in Refs. 
MM- For general fluids, we define quantities such as 


if 

= v° + vP, 

(20) 

p c 

= pA + 8), 

(21) 

p c 

= p + 5p. 

(22) 


They are respectively the velocity, density and pressure of 
the collapsed region. The background velocity expansion 
is given by and is associated with the Hubble’s law. 
Peculiar motions are denoted by iP. The total density 
within this spherical region under collapse p c is written 
as the sum of the background density and the overdensity 
fraction 5p. The same happens to the pressure definition. 

The rate at which the overdense region expands reads 

h = H+ 6 -(l + z), (23) 


Using the model parameters range which will be adopted 
in this work (10~ 6 meV< m x < 10 4 meV;10 -12 fm< 
l s < 10 12 frn) we calculate that 

k«2x (1CT 43 - 10 27 ) eV cm 3 , (18) 

which indicates that the Thomas-Fermi approximation 
can be adopted. 

Another comment about the justification of the use of 
the Thomas-Fermi approximation relies on the fact that 
we are going to focus on the largest cosmological scales. 
For example, in the Fourier space density perturbations 
are affected by the quantum pressure contribution pro¬ 
portionally to fc 4 while usual pressure contributions mod¬ 
ifies the evolution of the density contrast (which will be 
defined soon) according to k 2 [27|. Therefore, the quan¬ 
tum pressure corrections could be relevant for the very 
small scales. Besides, in the top-hat spherical collapse 
the density of all fluids inside the spherical overdense re¬ 
gion is homogeneous {39] and the anisotropic pressure 
m should be zero. 


where 0 = V • v?. 

Energy conservation is also required for the collapsing 
region. Therefore, each component i obeys a separate 
equation of the type 

Si = —3H(c 2 ff. — Wi)Si — [l + Wi + (1 + c e//i)^i] 

(24) 

where the energy density contrast is defined as 



and the effective speed of sound is computed following 
c 2 ff. = ( Sp/Sp)i. Note that over dot means derivative 
with respect to the cosmic time t. 

The dynamical evolution of the homogeneous spherical 
region will be governed by the Raychaudhuri equation 

q2 

9 + H9 H-= —47rGa (8 pi + 3 Spi) . (26) 

‘An ^' 
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For a cosmological model composed by N distinct flu¬ 
ids one has to solve N + 1 equations. One of the type 
l24l for each fluid and, since we adopt the top-hat profile, 
one single equation for the velocity potential 9 which is 
sourced by the density fluctuations of the N fluids. 

Since we will use the standard ACDM universe as our 
reference model here we show its equations for the spher¬ 
ical collapse. Both the baryonic and the dark matter 
component are assumed to be pressureless fluids. There¬ 
fore, we can write down 


and again, using the continuity equation we find 


Px — 


Pert 


(1 + W C rt)( ) 3 — W C rt 


Z — Zcrt ? 


(33) 


where p x is a continuous function at z cr t and w cr t = 
-yZL = cr 2 . At this point the continuity of the pressure 
(see discussion in 0 ) sets 


CT Pert — ^0 Pert ^ Pert — i (34) 

U 0 


5 b = -(l + S b )~, (27) 

a 

8dm — — (1 + 8dm) (l + O’ 2 ) —, (28) 

' a 

q2 

9 + H9 + — = — 4irGa [ pb8k + Pdm8dm{ 1 + o - ”)] . (29) 

Note that there is no equation of clustering of the cos¬ 
mological constant since it is treated as a background 
quantity. Therefore, it influences this set of equations 
only via the expansion rate H = H(ph, pd m , A). In or¬ 
der to numerically solve (I27H29D one usually specifies the 
initial conditions for 8b, 8dm and 9 at the redshift of de¬ 
coupling Zdec ~ 1000 from which one can treat baryons 
as an independent fluid. 


IV. ABRUPT PHASE TRANSITION 


The temperature T crt sets the beginning of the BEC 
phase transition. This is in fact a process which takes 
some finite time At until all the normal dark matter has 
been converted into the BEC phase. As estimated in 01 
At is of order of 10 6 years. Although the latter value is 
parameter dependent, it is in general indeed an almost 
negligible fraction of the universe’s lifetime. Therefore, 
the assumption that at z cr t there is an instantaneous con¬ 
version to the BEC phase seems to be plausible and it 
will be considered in this section. 

For z > z C rt the dark matter equation of state calcu¬ 
lated in © reads 

Pdm — ® Pdm , (30) 

where cr 2 = (u 2 ) /3c 2 . Applying this to the continuity 
equation one finds 

( l + z \ 3 <"- 2 ) 

Px = Pert ( i I ) ! Z ci z crt , (31) 

\ 1 + z crt J 


where z cr t is the redshift at the transition point and 
Pert = p(z C rt)- 

For z < z C rt the effective equation of state of the BEC 
dark matter is 


27 tH 2 I s 
m x 


(32) 


which, of course, depends on the model parameters. 
From this definition, 


_ Pert 

X (1 + Ldcrt){ 1 + Z C rt ) 3 ~ Wert 


(35) 


resulting in 


(1 + Zcrt) 3 


£lcrt 

fl x o 


’ ^ert 


1 + Wert 


(36) 


where Q x o = 0.25 is the today’s fractionary dark matter 
energy density parameter. The critical temperature at 
the point of Bose-Einstein condensation is 


Tcrt — 


2t tH 2 p 2/3 


ert 


<(3/2 ) 2 / 3 m x /3 k B 


= 6.87 


X 

( m x y /3 ( 

V1 meV J \ 


27 ih 2 


41 1/3 


a 


LC(3/2) 2 fc| ll 

2/3 


(37) 


3 x 10 


-6 


1 ftm 


- 2/3 


eV, 


where <((3/2) is the Riemann zeta function and Ub is the 
Boltzmann constant. 

Note that before the phase transition we have c 2 = 
cr 2 = w C rt- After this point, the equation of state pa¬ 
rameter and the adiabatic (c 2 = dp/dp) speed of sound 
associated to this fluid reads, respectively, 


w x (z) = u 0 p x {z ), c 2 3x = 2u 0 p x {z) = 2 w x (z). (38) 

Concerning the perturbed region the effective speed of 
sound is actually given by the expression 


„2 

C effx 


8p x 

8Px 


P° x -Px _ (l + aj 2 -! 
P c x-Px X 8 X 


w x (2+8 x ), 


( 39 ) 

from which one can expand for small values of 8 finding 
c 2 jj —>■ c 2 as expected. 

Since a crucial issue in this model is the determination 
of the moment at which the transition happens in Fig. |T| 
we show the dependence of z cr t on the model parameters 
m x and l s . This figure is numerically done after solving 
the equality proposed in (15C>1) . A giving z cr t value repre¬ 
sents a curve in the m x versus l s plane. The solid line sets 
the parameter values for which the transition happens to¬ 
day at z = 0. Therefore, only for the parameters values 
below the solid line the BEC dark matter model is able 
to leave some imprint on the observations. Note, for ex¬ 
ample that the configuration ( m x ,l s ) = (10 -4 meV , 10 5 


Px = u o P 2 X . 


u 0 = 
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fm) is an acceptable one. However, is this case, it would 
be impossible to probe the bosonic nature of dark matter 
since the transition will happen in a far future. On the 
other hand, over the long-dashed line the transition hap¬ 
pens at the time of photon-baryon decoupling. In prin¬ 
ciple, z cr t < 1000 is also allowed but its possible effect 
on the primordial CMB anisotropies is still not clearly 
known. Although this issue has not yet been investigated 
in detail we keep for convenience 0 < z cr t < 1000 where 
we can consider a matter dominated universe-apart from 
late A effects-and pressureless baryons. This redshift 
range corresponds to the gray region in this plot. The 
short-dashed line corresponds z cr t = 10 and it is shown 
to guide the reader on how z cr t evolves in this plane. 

We also show in this figure the usual range for axion 
masses 10~ 3 meV < 77i ax ion < ImeV. Taking typical axion 
scattering lengths < 10~ 16 fm, Fig. 1 estimates correctly 
that the axion condensation happens indeed very early 
in the universe history. In our work we are not advocat¬ 
ing in favor of any specific DM particle candidate. But 
in particular it is desired that most of the successes of 
the standard CDM paradigm should be kept. Indeed, it 
has been realised lon g ag o that axions are very promising 
candidates for CDM [43}. Therefore, the existence of such 
particles exemplifies the validity of our approach since it 
guarantees the non-relativistic behavior of the DM com¬ 
ponent before the phase transition takes place. Of course, 
there is no direct relation to actual CDM axion models, 
which condensate much earlier in the universe history, to 
our approach. Notice also that axions are characterized 
by an attractive self-interaction. We just use them as 
instance of CDM light particles. At the same time, our 
approach also relies on the fact that before the transition 
we are dealing with CDM like particles. Therefore, one 
should avoid to keep in mind the use of lighter particles 
since they would be associated to warm/hot dark matter 
models. 

The meaning of the mass of the dark matter particle 
is quite clear. But in the cosmological context what does 
the scattering length l s mean? 

Typical BEC experiments work with values in the 
range 10 6 fm< l s < 10 9 fm. For these values the con¬ 
dition z cr t > 0, i.e., assuring that the transition has 
already occurred, is satisfied for masses m > 10 meV 
and 100 meV, respectively. Of course, usual BEC experi¬ 
ments with atoms cannot guide us in our search for viable 
DM parameters. However, we also note that by extrap¬ 
olating the contours to lighter particles, as for example 
ultra-light masses of order m ~ 10~ 22 meV, z cr t > 0 
requires almost negligible l s values which can be much 
smaller than the Planck length ( l p ik ~ 10~ 2 ° fm). 

It is also worth noting that the space parameter indi¬ 
cated by the gray region is consistent with the stability of 
BEC dark matter halos as calculated in Ref. [-44}. How¬ 
ever, see also a related discussion on the non-stability of 
BEC halos in Ref. [45[. 

In order to solve for the evolution of the perturbed 
quantities during the collapse we adopt the following 



FIG. 1: The redshift of the phase transition ( z cr t ) in the 
parameters plane l s x m. The solid line sets the parameters 
in which the transition happens today at z = 0. The long 
dashed line sets the parameters for which the transition takes 
place around the decoupling time z cr = 1000. The axion mass 
range is shown only for the sake of comparison. 


strategy. We solve numerically the ACDM equations tak¬ 
ing initial conditions at a redshift Zi = 1000 and with the 
values Sdm(zi) = 3.5 x 10 -3 , Sb(zi) = 10 -5 and 9(zi) = 0 
EH S3] . These values represent the standard amplitudes 
in the linear perturbation spectrum associated to today’s 
clusters scales around the decoupling time. Indeed, the 
top-hat profile remains appropriate for such scales. No¬ 
tice that clusters scales collapsed at low redshifts and 
therefore already within the BEC dark matter epoch. 
Smaller scales which have collapsed before the BEC phase 
transition will preserve the CDM structure and only dif¬ 
ferences in the final virial configuration would exist which 
is not the scope of this work. With such initial conditions 
this set of equations is evolved until the critical redshift 
z C rt■ At this point, the quantities 6dm{z C rt), Sb(z cr t ) and 
6{z cr t ) are used as initial conditions for the BEC dark 
matter equations, which uses [39[ from the critical red¬ 
shift to z = 0. 

We have studied in great detail the parameter space 
m x and l s and although the BEC dark matter model in¬ 
deed yields to a distinct dynamics at nonlinear level, this 
difference is, in practice, almost negligible. We show in 
the left panel of Fig. [2] this feature where the expansion 
of the collapsed region is shown. The solid red line repre¬ 
sents the standard cosmology while the dashed black line 
was calculated for a mass m x = 20 meV and a scatter¬ 
ing length l s = 10 6 fm. With this choice the transition 
occurs at z cr t = 3.19 as seen in the vertical dashed line. 
Both curves are in practice indistinguishable. The effec- 
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FIG. 2: Expansion rate (left) and effective speed of sound 
(right) of the collapsed region. In both plots we have m x = 20 
meV and l a = 10 6 fm. 


The dynamics shown in this section was also developed 
for the first time in Ref. [Tsl |. 

As mentioned in the last section the estimated duration 
At = t(zBEc) ~ t(z cr t) of this transition is of order At ~ 
10 6 years, that is a small fraction of the universe’s lifetime 
tjj ~ 10 10 years pj|. However, At depends on the model 
parameters l s and m x . We calculate here again At for 
some values l s and m x and plot it in Fig. [3] In the right 
panel of this figure, there is a maximum value At max = 
3.4 x 10 9 years assuming, for instance, a mass m x = 1 
meV and l s ~ 3.1 x 10 2 fm. There are of course other 
combinations of l s and m x which produces similar At 
values. The lower values for At we have found are ~ 
10 6 years. Therefore, this analysis shows that contrary 
to previous estimations, the phase transition can last a 
non-negligible fraction of the universe’s lifetime. It is 
worth noting that recently Ref. [46| has pointed out the 
preferred values m x ~ 10 -3 meV and l s ~ 10“ 7 fm which 
according to our Fig. 3 maximize the duration of the 
phase transition. 

As we will see below, the background dynamics and 
the evolution of the perturbation for the smooth phase 
transition differs significantly from the abrupt case stud¬ 
ied in the last section. Then, one can expect that now 
we can observe some distinguishable feature of the BEC 
dark matter nonlinear collapse. 

Let us now develop the dynamics during the smooth 
phase transition. Before the transition starts, we have 
the same dynamics of a isotropic non-relativistic gas, as 
described in the last section by equations m and ED- 

During the phase transition we can define the fraction 
of converted dark matter as 

m= p(z) ~ pcrt , (40) 

PBEC — Pert 

where p(z) is the dark matter density along the transi¬ 
tion, p crt is the dark matter density before the transition 
and pbec its value afterwards. The function f(z) is de¬ 
fined in such a way that at z cr t we have f(z cr t) = 0. 
When the dark matter has fully converted to the BEC 
phase / (~bec) = 1- 

Using (14011 into the continuity equation and integrating 
it from z cr t to z > zbec we find 


five speed of sound is plotted in the right panel of Fig. 
[2j This shows the reason there are no significant changes 
in the evolution. We remark again that this result is not 
due to the specific choice m x = 20 meV and l s = 10 6 fm. 
It is a general feature of the model. 


V. SMOOTH PHASE TRANSITION 

We deal now with the situation in which there is a 
gradual conversion of “normal” dark matter into the con¬ 
densed phase which starts at a redshift z cr t and is finished 
at a redshift zbec- This is indeed the more realistic case. 


“ert L J 

Then, the dark matter density evolution becomes 


IT z 

p x = p ^{TT^r t 


3(1+<t 2 ) 


Z A Z C rt '■ 


P\ ~ Pert ^ 1 + (1 + U cr t) 
Zcrt > Z > zbec ; 

Px — Po- 


1 + z 


-1 


(42) 

(43) 


(1 + z ) 3 


(1 + W 0 ) — Wo(l + z) 3 


1 + Z cr t J 
, Z < Z B EC • (44) 



















0.4 



FIG. 3: The phase transition time length At as a function of 
the models parameters l s and m, where we fixed a 2 = 3x 10 -6 . 


We still have to the determine the redshift zbec 
when the phase transition is over. With the condition 
/(zbec) = 1 inserted in (TUI) we find 


^BEC 
H crt 


1 + (1 + tO cr t) 


1 + Z cr t 
1 + ZBEC 


3 

= (1 + U} cr t) , 


(45) 

and using 0 = zbec in the expression (l44l) for the con¬ 
densed dark matter density we have 


^bec _ n° rt (1 + z bec) 3 
fieri (1 + Wo) — Wo(l + ^BEc) 


(46) 


Eqs. © and CCH) can now be solved, leading to a solu¬ 


tion for zbec and S4bec- 

As said before, during the phase transition both non- 
condensed and condensed dark matter coexist and the 
dark matter pressure is constant having the same value 
for both components in the interval zbec < z < z CT t, as 
given by Eq. [34] We will assume that the same happens 
for the collapsed pressure p c . This allows us to find the 
constraint 


1 + 6* Tt = (1 + ^ rt ) 2 , (47) 

where we used the expression p x (z) = p a {z) + Pb(z ), 
which compared with Eq. m allows us to identify 
p a {z) = Pert(1 — /(z)) as the non-condensed dark matter 
density and ps(z) = pbec/(~) as density of the con¬ 
densed state. 

The continuity of dark matter fluid pressure enable 
us to treat both components as one single fluid also at 
perturbed level. In this case, the effective fluid sound 
velocity during the phase transition becomes 


c 


2 

e ffx 


Pert ~ Pert 
Px^x 


_ p p CTt 


Px 


<5crt 

S x 



(48) 


where w(z) = w cr tPcrt/ p x (z) is the equation of state pa¬ 
rameter for the dark matter fluid during the phase tran¬ 
sition. After the phase transition is completed, i.e., when 
z < zbec, the effective fluid sound velocity of the BEC 
dark matter will be 


.2 

C e//x 


Pert Pert & Pcrt^crt 


Px^x 


Px^x 


= w(z)(2 + 5 x ), (49) 


where ui(z) = w crt p x (z)/p cr t is the equation of state pa¬ 
rameter for the dark matter after the phase transition. 

Since the velocity dispersion cr 2 for the dark matter 
particles before the BEC phase transition is small the 
same assumptions made on z cr t in the previous section is 
still valid here, and we will consider values for the model 
parameters ( m x ,l s ) such that 0 < z crt < 1000. We will 
also consider only cases where zbec > 0. 

This set of equations is evolved until the critical red- 
shift z C rt assuming the same initial conditions at z, as be¬ 
fore. At this point, the quantities 1 5dm.(z cr t), 8b(Zcrt ) and 
9{z C rt ) are used as initial conditions for the phase tran¬ 
sition perturbed Eqs. © and © with the suitable 
background parameters. This set of equations is again 
evolved until the zbeCj and the quantities <Wi(2bec)> 
^(zbec) and 0(zbec) are used as initial conditions for 
the BEC dark matter perturbed equations. 

In Fig. 0]we show the expansion of the collapsed region 
for the smooth phase transition model, where the solid 
red curve represents the standard ACDM model and the 
black dashed curve represents the BEC model, for m x = 
20 meV in the left panel, m x = 10 meV in the right 
panel and l s = 10 6 fm in both cases. The dashed vertical 
lines show the initial and the final points of the transition 
phase. In the left panel, z cr t = 3.19 and zbec = 1-43 and 
Zcrt = 1.10 and Zbec = 0.45 in the right panel. These 
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Z 



Z 

FIG. 4: Expansion rate of the collapsed region for the smooth 
phase transition approach. 


intervals correspond to 2.40 x 10° years and 3.38 x 10 9 
years. As in the abrupt transition model there are no 
major difference between CDM and BEC dark matter. 

The evolution of the non-linear density perturbations 
are shown in Fig. [5J where 5dm = 5pdm/pdm is the dark 
matter density contrast. Again, the red curve represents 
the standard ACDM model, while the black dashed curve 
shows the behavior of the BEC model for m x = 20 meV 
in the left panel, m x = 10 meV in the right panel and 
4 = 10 6 fm in both cases. The curves are again indistin¬ 
guishable. 

The redshift of turnaround 2 ta is the one which marks 
the moment when the perturbed region starts to decrease 
its physical radius. This happens when h = 0, i.e., z ta = 
z(h = 0). For ACDM model z t A a CDM = 0.2113 and for 



Z 



Z 


FIG. 5: Dark matter density contrast for the smooth phase 
transition approach. 


the cases seen in both panels of Fig. [4] we have |z EEC — 

4 CDM l«io- 4 . 


VI. CONCLUSIONS 

We have studied the nonlinear clustering properties 
of the Bose-Einstein dark matter model. In this sce¬ 
nario, bosonic dark matter particles are able to undergo a 
phase transition as their temperature reaches the critical 
one T crt which corresponds to some critical redshift z cr t- 
The main questions here are: i) how does z cr t depend 
on the fundamental model parameters m x (the particle 
mass) and l s (the scattering length)? and ii) what is the 
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background and perturbative dynamics during the phase 
transition? 

Fig. [T] shows in detail the expected degeneracy of z cr t 
values in the l s x m x plane, i.e., for a given z cr t , there 
are many admissible parameter configurations. This re¬ 
sult identifies the parameters values for which z cr t > 0 
and therefore are able to leave imprints on large scale 
structure observations. At the same time, if the actual 
parameters values of the BEC model lie in the region 
z cr t < 0 then the bosonic nature of the dark matter par¬ 
ticles cannot be accessed via cosmological observables. If 
the present model is employed for BEC phase-transitions, 
ultra-light candidates (m x < 10 -22 eV) would only lead 
to possible observational imprints for l s of order of the 
Planck length or smaller. Note that this claim is limited 
to the fluid description used here. Recent calculations on 
the full dynamics of the ultra-light axion scalar field sho 
that there are indeed possible observable imprints in the 
cosmological data 01 • 

Our strategy was to identify specific signatures of the 
BEC dark matter nonlinear clustering. Since there is 
a positive pressure associated to the BEC dark fluid 
one can expect that the corresponding effective speed of 
sound will modifies somehow the agglomeration rate. We 
tried to understand this process via both the abrupt and 
the smooth phase transition approaches. In the former 
scenario the dark matter dynamics changes suddenly at 
z cr t■ In the latter, there is a continuous conversion from 
the “normal” to the BEC phase. Although we showed 
that the smooth transition can indeed last quite a signif¬ 
icant fraction of the universe lifetime. Then, it seems 
that this case could, eventually, lead to a remarkable 
dynamics. However, in both approaches of the phase 
transition we could not identify any relevant difference 


between the BEC model and the standard CDM model. 
This is mostly because the model parameters leading to 
z C rt < 0 produce almost negligible c e // values. On one 
hand, this guarantees that the nonlinear clustering pat¬ 
terns of the BEC model at large scales are very similar 
to the CDM model. We have provided a theoretical con¬ 
firmation for the recent numerical results of Schive et al 
[33j which claims that the differences between BEC DM 
and standard CDM appears only in the internal struc¬ 
ture of DM halos rather than on the cosmological large 
scale distribution. On the other hand, this eliminates 
the cosmological nonlinear perturbative study as a possi¬ 
ble technique to probe the bosonic nature of dark matter 
particles. It is also worth noting that the typical value 
for the critical overdensity for collapse S c = 1.686 re¬ 
mains unchanged for the BEC parameter space probed 
here. Perharps, this conclusion is in part due to the fact 
we have assumed a simple version of the first order phase 
transition of the BEC DM model. Taking properly into 
account, for example, the latent heat released during the 
transition and the resulting dynamics associated to the 
nucleation of the new bubbles we could end up with a 
very drastic effect on the non-linear clustering. We will 
leave this analysis for a future work. 
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